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ABSTRACT 

A variational principle and a finite element discretization technique were used to 
derive the dynamic equations for a high speed rotating flexible beam-mass system embedded 
with piezo-electric materials. The dynamic equation thus obtained allows the development 
of finite element models which accommodate both or the original structural element and 
the piezoelectric element. The solutions of finite element models provide system dynamics 
needed to design a sensing system. The characterization of gyroscopic effect and damping 
capacity of smart rotating devices are addressed. Several simulation examples are presented 
to validate the analytical solution. 

INTRODUCTION 

Structural monitoring and control using smart sensing materials has become of 
importance in recent years due to the rapid development of large flexible structures and 
flexible mechanical systems. These materials and structures have their own sensing, 
actuating, tuning, controlling and computational abilities (Gandhi and Thomas, 1989). 
Typically, smart materials and structures are distributed- para meter systems operating 
under a variety of service conditions and having a theoretically infinite number of vibration 
modes. Current design practice is to model the system with a finite number of modes and 
to design a sensing system using lumped parameter approach. Truncating" the model may 
lead to performance trade-ofT when designing a control system for distributed parameter 
systems. Lumped parameter approach is generally acceptable for sensing applications due 
to its nature of simplicity nnd case of implementation. Although significant progress has 
been made in the recent past in the development or smart materials and structures 
featuring piezo-electric materials (Bailey and Hubbard, 1985; Plump et.al, 1987, Crawley, 
19887), shape memory alloys (Miwn, 1985; Yacgcr, 1984; Rogers and Robcrtshaw, 1988), 
electro-rheological fluids (Gandhi et al, 1989; Choi ct al, 1989; Gandhi ct nl, 1989), and 
optical fibers (Morikawa, 1985; Rogowski, 1988), very few work hove been done to 
characterize the dynamic behaviors of these devices, especially when H is used for high- 
speed applications. 
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Emphasis is plnccd on the accurate modeling and characterization of structural 
parameters of sensing devices for flexible structures. The analysis of a beam attached to 
a rotating base is a subject of interest to many researchers because numerous structural 
configurations such ns spinning satellites (Laurcusnn, 1976; Kane et al., 1987) and flexible 
robots (Cannon and Schmitz, 1984; Mitchell and Bruch, 1988; Yang and Donath, 1988) fail 
into this category. The analysis of these rotating dynamic structures with payloads is quite 
different from those of statinnaiy structures due to the inertia of gyroscopic effect at high 
rotating speeds. 

Piezoelectric materials are media which develop mechanical strain when subjected 
to an electrical field, or conversely, they develop an electrical field when subjected to 
mechanical deformation. Their inherent high power-to-weight ratio makes them ideal 
candidates for embedding piezoelectric materials in traditional structures for vibration 
sensing and control. Crawley and de Luis (1987) studied the effect of a beam with bonded 
piezoelectric sensors. Plump et al. (1987) used a piezoelectric film to enhance the damping 
ratio of a cantilever beam. Tzou applied a piezoelectric film os an active vibration in a 
flexible structure (1987) and an active vibration isolator and exciter (1989). This paper 
focuses on the development of analytical models for dynamic characterization of a high- 
speed rotating flexible beam-mass system with embedded sensing system. The design of 
such a sensing system featuring piezoelectric materials is addressed. Several simulation 
examples are presented to validate the analytical solutions. 

DESIGN AND ANALYSIS 


Basic Assumptions and Coordinate Systems 


To derive a simple yet effective sensing model Tor the physical system of interest, 
several assumptions are imposed here, namely: 

(1) Large Overall Rigid-Body Displacements with Small Elastic Deformations: This 

assumption is valid for a mechanical system rotating at speeds less than one 
thousand revolutions per minute with low payload. 

(2) Negligible Gravity Effect: For high-speed rotating space structure, the gravity effect can 

Ik? ignored. 

(3) Negligible Geometric Stiffening Effect: When rotating in plane, the magnitude of the 

axial displacement Is much smaller than that of the transverse displacement. The 
effect of geometric stiffening is negligible. 

(4) Plane Stress Condition: Since the beam is thin, the stress variation through the 

thickness is negligible. 

(5) Average Material Properties: The average material properties of smart beams is used. 

Since the piezoelectric film is relatively thin, the isotropic nluniinum beam plays a 
dominant role in contributing to the overall beam deflection. For structural 
monitoring, this assumption is valid. 


_ Mi + pA 

P, t - 

E _ e a : 

** f 


(i) 


where A t is the cross sectional area of AL layer, A 2 is the cross sectional area of PVF, layer, 
i, is the moment or inertia of Aluminum layer about z axis, 2, and L is the 


l*o coordinate systems, one being the global X-Y and the other local x-y, are 
inlrodnced to describe the dynamic system. The local coordinate system moves with the 
rigid body configuration of the link. Figures 1 and 2 show the schematic of the deflected 
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Figure 1: Schematic of a 
deflected beam element in 
local and global coordinate 
axes. 


Figure 2: Schematic of 

rotating beam-mass system 
showing local and global 
coordinate axes 


beam in both local and global coordinates. In classical small dellection beam theory, the 
displacements of an arbitrary point P(U, W) on a beam can be expressed in terms of 
homogenous coordinates as 


[Hi [10 0 
1/ - 0 -y 1 

[ij 0 0 0 



( 2 ) 


where x and y are the local coordinates of point P in the undeformed state. U and W are 
the local coordinates of point P after deformation, u, w, and w* are the axial, transverse and 
tangential displacements of point P, respectively, {g} is the displacement vector of point 
P in homogenous coordinates. The displacement vector, {g}, can be expressed in terms of 
nodal displacement vector {p} ns {g} = N <p>, where {p} = [{q} T | I ] T , {q} = [w, w 4 ' u, 
Wj w; iij) T , and N, the shnpe function matrix, is given as 

N, N 2 0 W, JV 4 0 0 
N[ 0 A/J ^ 0 0 

N - 

0 0 N s 0 0 0 

0 0 0 0 0 0 1 




where the shnpe functions, N,, N 2 , N 3 , N 4 , are the Ilcrmile polynomials used for a beam 
element. N $ and N 6 arc the shape funclions for a Imr element in nxinl loading [24]. Note 
that N,‘ = dN/dx. These shape functions are reported as 


N, - l-3(i) 1 .2(i) 1 ; W 2 - 

n, - N > - - <■{) 


Substitute the expression of {g} into equation (2) to get 
[W U \] T -BN[p\ ( 5 ) 

and describe point P in terms of the global coordinates ns 
\X Y \} T -RBN[p\ ( 6 ) 

where R is the transformation matrix that relates the local and global coordinate systems. 
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The.expressJon for R is 
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,he ,wca, coordlnat * a *«* «nd 0 Is the angle 
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Elemental Equation of Motion 

equation ore used to derive' foe ^ment*' for 0 d yn«nilc system and the Lagrange 
formulation is to discreti^^ fP «■ ‘-e finite eK 

of foe Lagrangian of its constituting element, i \ he gronglan °ra system is the sum 
the Lagrange equation can be presented as ' 1 ° f kine ‘ ,C and P 0,e,,tin, energy, 
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-» V s? "r - * *. 

Of the rotating beam and the tip mass attached at the'cnri”" T " b f in ,* ,he klnetic energies 
energy of the beam element, {q} and (Q V are , he . ’ U ,s ,he potential 

respectively. The expressions of M, hI M M KM* f?* COOrdlnaU * «nd forces, 
The detailed derivation of these matrices are%£ “dln^uVft 8 ' Ven be,<w - 
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where I and A represent the moment of l ™ tih ^ JlveV m bW,m ' 
respectively. The noniero terms of M lni , M, m , M_ » na V,« 

IMJm - “J> - 

t"Ju - - <r.«> tow ' ' 

t«Ju * "J 1 ' <*.«**' * '.* ,w ^ 1 

(MJ U - WjW}<z*J ♦ - 1) 

1V,J„ - M.(4 ♦ >0«j. - ’■^ - o. * J > xn ^ (15) 
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Note that the system mass matrix is the sum of the mass matrices of the tip mass 
and the beam element, l.e., M 2 - M^+M^, and so on. 

System Equations and FEM Solution 

slate, the elemental equation, of mollnn are t™eraled. '■™^ emb| , J 

»l“‘"r I, oblalnml b, .be In.e, union at.he 

v oUage , ore calculated using the result of finite element solutions. 

Stress-Strain Relationship 

The longitudinal strain €„ Is given as 


, t . 21 . * - rS - i« -y » u 






(17) 


Stress at the same point can be obtained by using 
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Dynamic Monitoring System 

From these equations, strain and stress at any arbitrary point, p(*» y)> on *k*J*f m 
element are computed. The voltage, V(x,t), generated by the piezo strain is calculated by 

, v (18) 

K(jc,0 - e f (x,t).(hjd n ) 

where E Is the modulus or elasticity, h is the thickness or the layers, V is the strain induc^d 
voltage, h, is the thickness »r the PVF, layer, and d„ is the appropriate piezoelectric 

constant, and 

*.*. + £ A fjH (19) 

*' {x) — * 

The sensor configuration used in this study is a layer of PVF,, polyvinylldene 
fluoride, bonded to one side or the rotating beam. Fig. 3 shows the cross sectional view of 
the beam with PVF, layer. PVF, Is a polymer that can be polarized or made piezo- 
electrically active through appropriate processing during manufacture (Bailey an u a , 

Si ll nonpolarized Elm, PVF, i. » 

fnriii PVF is essentially a tough, ncxiblc piezoelectric crystal. Polarized PVF, is 
commercially available as thin polymeric film having a layer «r nickel or » un,lnum 
deposited on each lace to conduct a voltage or applied across its races In y-direction which 
results in a longitudinal strain in x-directlon. This is the d„ component or the piezoelectri 
activity, ir PVF, is polarized biaxially that would strain in both the x and the z directions. 
For this study we consider uniaxial PVF, only. 

Fig. 4 shows the smart beam configuration. The longitudinal strain, Is obt ^ i “f d 

by Torce equilibrium in axial direction by finite element formulation. The strain, « p , in the 
J f ia tUm. f>VIT This mmnined 


effect of gyroscopic, coriolis and other inertia are recorded by the piezo voltage. This 
voltage con be used as an monitoring index. When the Index v»lue exceeds the imposed 
constraints or material’s natural constraints, the control action, i.e. the direct piezoelectric 
effect, can be triggered to take place. 

Control Algorithm 

The piezoelectric strain creates the net force in each layer acting as the moment arm 
from the midplane ot the layer to the neutral axis of the beam, producing a bending 

moment (20) 

TUA - £,*,>*1 K*A>*V B 1 

where b is the width of the beam, and D is the location of the neutral axis of the composite 
beam given by 

( 21 ) 
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Performing some algebraic manipulations to yield 


( 22 ) 
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where c is constant for a given beam material and 
per volt. If the material property and geometry 
length, c is a function of x. 


geometry expressing the bending moment 
of the composite beam change along its 


Combining ,h. .born »»b * 

equations of motion for transverse vibrations, w(x,t), of the compt 

equation becomes 


Ay 1 


. &W 

,A— 


0 ; 


for 0<x<L 


(23) 


and the boundary conditions are 
B*\ - M,^l t * c.^ 

&> a* 1 ** 


(24) 



1 #w 3 w (25) 

of the first variation of the functional gives 
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where -2— | t is the angular velocity at the tip of the beam. The control voltage Is chosen 

to generate a bending moment that opposes the angular motion at the tip of the beam. The 
geometry oMhe piezoelectric layer can be tailored to obtain the necessary control function. 


NUMERICAL EXPERIMENTS 


Several numerical simulations are conducted to illustrate the effectiveness of the full 
beam formulation of the beam-mass system. A beam-mass system made of rectangular 
aluminum materials with an embedded PVF, layer Is used for the study. The dimensions 
and material properties of the beam are given in Table 1. The beam is subjected to a spin- 
up maneuver (Kane, 1985) prescribed by 


; omT 

M'-fl : ' > T 


(29) 


where 6 0 is the steady stale angular velocity. In order to characterize the dynamic 

behavior of the system at different rotating speeds, we set the steady state speeds to 300, 500 
and 700 rpm, and time constant, T, is 1 seconds. The simulation is carried out for a period 
of 5 seconds. 

Table 1 - Geometric parameter! and material properties of U>c wart be— 



Aluminum 

PVF, 

Corny— iW 
Iks— 

LtngUi (L) 

20 

20 

20 

Thkkam (H In) 

.9999 

.0011 

( 

WMlh (w; In) 

.IS 

.25 

.25 

Ifcndtr (ptf » tW) 

145 x 19 4 

1.4* X 10 4 

145 x I0 4 

Yaangt'i M«4wl» (K| H> 

l I I0 7 

tf i l» 5 

I x to 7 

SUIk Ptc— -rkclrfc 


!M s ir 10 



Figures 5 and 6 show the transverse displacement, w, of the piezo-aluminum beam- 
mass system rotating at three different speeds. The transverse tip displacement of the beam 
is computed and recorded for both formulation with and without payload mass inertia. 
Both of the maximum transient tip displacement and the steady-state tip deflection are 
proportional to the size of the mass attachment. The contribution of gyroscopic terms is 
computed and recorded in Table 2. A maximum 4.35% contribution of gyroscopic Inertia 
is observed at a speed of 700 rpm. The contribution becomes significant when the mass 
attachment is increased. The contribution of gyroscopic Inertia becomes more pronounced 
when the speed is increased. 

When 10% of a equivalent beam mass Is attached to the tip of the rotating beam, the 
tip strain is increased at least two times. Figures 7 and 8 show the longitudinal tip strain 
at three different speeds. The corresponding longitudinal stresses are presented in Figures 
9 and 10. These plots are similar to those of Figures 7 and 8 with an amplification in 
magnitude by the factor of an equivalent modulus of elasticity. The tip strain induced 
voltage of the piezoelectric layer at three speeds are presented in Figures 11 and 12. It is 
observed that the piezo voltage introduced by tip strain of no payload Is In the range of 50 
to 170 volts. The piezo voltage induced by the tip strain with a 10% tip mass attachment 
is in the range of 120 to 530 volts. The contribution of tip masses and speeds are 
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experimented to allow control voltage to be adjusted to the design range. The geometric 
parameters of both aluminum beams and piezoelectric layers are design variables to be 
tailored to satisfy the design need. 


Table 2 - Contribution of gjrrtweoplc term* at JilTcrtnt apetda (*) 



Tip Mau 

Angular Velocity | 

Attached 
lo Beam 

300 rpm 

900 rpm 

700 rpm 

Aluminum 

bummui 

•yetem 

0% 

0.92% 

1.02% 

2.06% 

10% 

0.67% 

1.47% 

3.20% 

I Plezo-alumlnum 

H beam -mail 

| ayiuni 

0% 

0.69% 

1.36% 

2.60% 

10% 

0.90% 

2.06% 

4.39% 


(*) percentage difference in eteedy-itate eotullon at drtfecant epeeda foTSf^and M 1 



500 sad 700 rpm 



700 rpm 


CONCLUSION 


A systematic finite element based design method is presented in the paper. The method 
allows the a high-speed rotating dynamic structure with embedded piezoelectric films to be 
designed and monitored. The gyroscopic effect introduced by different paylond at various 
speeds can be taken into consideration in the early design. The result of numerical 
simulations indicates that current appronch can be used Tor application in the sensing and 
monitoring of high-speed spinning space structures and flexible mechanical systems. The 
finite clement based method Is simple anil systematic. The dynamic characteristics of high 
speed rotating sli octurcs and ninchiiicry can be observed and used in distributed parameter 
models for control of such systems. 
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Figure 11 - Piezo voltage induced 
by the tip strain or plezo- 
aluminum beam rotating at 300, 
500, and 700 rpm 



by the tip strain of the plezo- 
aluminum beam-mass (10%) 
system rotating at 300, 500, and 
700 rpm 
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